Homogeneous Fermion Superfluid with Unequal Spin Populations 
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For decades, the conventional view is that an s-wave BCS superfluid can not support uniform spin 
polarization due to a gap A in the quasiparticle excitation spectrum. We show that this is an artifact 
of the dismissal of quasiparticle interactions V qp in the conventional approach at the outset. Such 
interactions can cause triplet fluctuations in the ground state and hence non-zero spin polarization 
at "magnetic field" h < A. The resulting ground state is a pairing state of quasiparticles on the 
''BCS vacuum". For sufficiently large V qp , the spin polarization of at unitarity has the simple form 
Our study is motivated by the recent experiments at Rice|l[ which found evidence of a 
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homogenous superfluid state with uniform spin polarization. 



What happens to an s-wave BCS superfluid when the 
number of up and down spins become unequal? For 50 
years since the invention of BCS theory of supercon- 
ductivity, this fundamental question remains unsettled. 
Recent experiments on degenerate Fermi gases of Li 
at MIT 2] and at Rice University y|, however, have shed 
lights on this long standing problem. 

According to mean field theory, the BCS ground state 
is highly resistant to chemical potential difference <5/i = 
2h between the two spin populations. The system can 
not develop spin asymmetry (i.e. with zero susceptibil- 
ity) unless h is of the order the energy gap A. It is 
found that for h > A/V2, the BCS state is unstable 
against the spin polarized normal stateQ. Exactly how 
the BCS state evolves as h increases, or more relevant to 
current experiments, as spin polarization increases, is a 
subject of controversy. There have been many propos- 
als. The most famous one is the so-called FFLO (Flude, 
Farrel, Larkin and Ovchinikov) state[j|; which is a pair- 
ing state with an oscillatory gap along a specific spatial 
direction. Recent studies of strongly interacting atomic 
Fermi gases 4], however, show that the FFLO state can 
exist only in an insignificant range of spin polarization, 
and that any amount of polarization will cause the gas 
to phase separate into regions of unpolarized BCS state 
and spin polarized normal state. In the following, to 
make contact between atomic Fermi gases and electrons 
in metals, we shall refer to the difference in spin densities 
(?i| — ri|) of a Fermi gas simply as "magentization" m, 
and h simply as "magnetic field" . 

In both Ref.dQ, the density profiles of different spin 
populations of 6 Li Fermi gas in the strongly interact- 
ing regime have been measured. The MIT data are for 
P > 0.1 (P = (n-f — nj.)/(n-| +nj,) = m/n) taken after the 
gas is released from the trap. The deduction of original 
spin density requires the understanding of the expansion 
dynamics of a strongly interacting Fermi gas - an intrigu- 
ing problem yet to be worked out. The measurements of 
the Rice group were performed in situ, including polar- 
ization < P < 0.1 (with accuracy 0.03). Although 
there are no explicit determination of temperature, the 
experimental conditions are similar to those in earlier ex- 



periments in the superfluid phase. It is then reasonable to 
think that the system is a superfluid. While both groups 
have observed similar spin density profiles at larger po- 
larization, the Rice group also has found evidence of uni- 
form spin density when spin polarization P < 0.1. This 
is a highly surprising because such a state is forbidden 
within mean field BCS treatment. (Although a recent 
Monte Carlo calculation suggests some features of such a 
state, the evidence is very indirect that it is hard to make 
a strong casepj). While more experiments are needed 
to further confirm this finding, it raises the fundamen- 
tal question of whether there is any mechanism at all to 
allow Fermion superfluids to accommodate uniform mag- 
netization. 

In this paper, we show that quasi-particle interactions 
(due to interactions left over when reducing the actual 
Hamiltonian to the BCS form) is a natural mechanism 
to generate triplet spin fluctuations. These fluctuations 
will lead to a non-zero response to "magnetic field" h, and 
hence a homogenous superfluid with uniform spin polar- 
ization. To realize this true ground state mathematically, 
however, requires taking the kind of conceptual steps for- 
ward as in BCS theory. It is well known that due to 
number conservation, BCS state cannot be obtained by 
performing perturbation theory on the normal state. It 
can only be realized by constructing a coherent state in 
number space to allow pair fluctuations. Likewise, due to 
the spin conservation, the true ground state in non-zero 
magnetic field can not be realized by performing pertur- 
bation of quasi-particle interactions on the BCS state. It 
can only be obtained by considering coherent states in 
spin space allowing spin fluctuations. 

A. The natural emergence of triplet excitations: 
Let us first ask a simple question: What would be the 
ground state of the BCS Hamiltonian if we insist that 
the system has a very small spin, (say S — 1), and 
with zero total momentum. The answer is to create a 
triplet pair with zero momentum on the BCS ground 
T | BCS), where 
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P k creates a coherent state of pairs in the number space, 
u k and v k are the coherence factors. Since \BCS) is 
annihilated by the quasi-particle operators a k = M k c k| + 
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and /3_ 



-w k c k | + w k c_ k |, we have 
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In other words, creating triplet excitations from the real 
vacuum is equivalent to creating quasi-particle pairs from 
the "BCS vacuum". Since attractive interactions cause 
fermion pairing, attractive interactions between quasi- 
particles will cause pairing between them. 

Bl. Interaction between Quasi-particles: The 
Hamiltonian of atomic gases is H = T + V, where 

T = Ek [(& - M c kT c kt + (& + /l ) c ki c k|]> Ck = e k - M) 
£k = h 2 k 2 /(2M), and V is pseudo-potential for two- 
particle scattering [||. V conserves both number and 
spin. In the BCS approach, one considers a ground 
state \BCS) (eq.JTJ) that allows fluctuation of singlet 
pairs. This leads to the reduced Hamiltonian H BCS — 



r+E k [A, 
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-/i.e.], where A is determined self con- 



sistently as A = 9 5(r)9 r [r(^ ; (R-r/2)V' T (R + r/2))] 0, 
<7 = 4Trh 2 a s /M, and a s is the s-wave scattering length. 
Explicitly, we have 
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where = y (ek A*) 2 + A 2 . We then have 



BCS — 



k 



[(£k - &)4a k + (E k + h)l3{p v 



E 



BCSi 



(4) 



where E BCS = X)k(£k — E^) is the energy of the BCS 
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state (eq.JTJ) 8]. The fact that eq.© and £ 
dependent of h means the \BCS) state is rigid against 
spin polarization, since m = —d(E BCS /VL)/dh = for 
h < A, where f2 is the volume of the system. 

To illustrate the essential physics of quasi-particle in- 
teraction, we consider the following model interaction 
V qp , (its origin will be discussed in section C): 

^ P -^E w(p,p')rpr p „ r p , = a P ,al p ,-/3_ p ,/3 p ,, 
p,p' 

(5) 

where W(p, p') is symmetric in p and p'. It therefore 
has the decomposition W(p, p') = Ei, (p^iUp') 
where io's are orthogonal, O -1 ^ u;„(p)u;*/ (p) oc 
Eq.© also implies that ^(PjP'j is odd in p, hence 
w i>(p) = — Wu{— p)- We shall assume W has at least 
one positive eigenvalue, i.e. —W has an attractive com- 
ponent. Finally, we shall assume that the eigenfunctions 
w decrease sufficiently fast at large wavevector so that 

h = E P K(P)I7£ P converges©. 

B2. A new ground state: Our Hamiltonian is now 
TL = Hbcs + V qp - Note that V qp conserves polarization 



M = N+ — N,. Since TL (for quasi-particles) has the 
same form as the pairing Hamiltonian for fermions, an 
attractive W will augment the "vaccum" (i.e. \BCS)) 
with quasi-particles pair fluctuations, which, as discussed 
in Section A, are triplet fluctuations. A non-zero h will 
then produce more |T than J, J, pairs and hence a non-zero 
polarization M. 

There are, however, major differences between pair- 
ing of fermions and pairing of quasiparticles. Since the 
spectra of fermions (£ k ) are gapless and those of quasi- 
particle (-E k ) are gapped, pairing has much more dra- 
matic effect on the former. Furthermore, in the absence 
of pairing, the ground state of fermions consists large 
number of particles (normal Fermi sea) whereas that of 
TL (i.e. BCS "vacuum" \BCS)) consists of no quasipar- 
ticles. Thus, the re-organization of normal Fermi sea 
caused by fermion pairing is much more dramatic than 
that of the \BCS) caused by quasiparticle pairing. Nev- 
ertheless, such reorganization is sufficient to cause the 
system a non-zero response to ft,, a property that does 
not exist in the BCS state. 

The analogy with BCS state also shows that the true 
ground state (with M at h < A) can not be obtained by 
performing perturbation of the spin conserving V qp on 
\BCS), which preserves the zero polarization of the BCS 
state despite h ^ 0. What is needed is to go to the grand 
canonical ensemble in spin space, allowing fluctuations 
in spin polarization, and replacing TL by the mean field 
Hamiltonian 
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BCS 
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with mean field 
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Spin symmetry is now broken. It is then possible 
for the system to produce magnetization to gain en- 
ergy from h. Eq.© can be diagonalized as TL = 

E P (4~ )A l A p + 4 +) K B p) + E G, with S round state 



\G) = II (u P + v P al p al) (u p + v p f3l p (3 P ) \BCS), 



where u k , v u , uu, vu are coherence factors 10]; Eq 
E P Up - [4 _) + 4 +) ]/2) is the energy of |G) and 
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£ (E p -[£<->+ 4 +) ]/2) <0. (9) 
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The operators A p = u p a p + v p aL p , B p = u p /3 p + v p f3_ p 
are the new quasi-particles with energies 



£ p a) = J(E p + ah^ + \r\ 2 , a = ±. (10) 
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Eq.© and 1|10|) show that spin fluctuations £ k increase 
the energy of the excitations but decrease the energy of 
the ground state. They also lead to a "magnetization" 

to = n- 1 Ek( c k T c k T - c -ki c -ki> = Ek(«k - ^k)> 



a{E^ + ah) 



(ii) 



B3. Triplet mean field and magnetization m: 

Evaluating (r ,) in eq.0, we have 



p',CT=± \ 



p' 



V(^p'+^) 2 + |C P 'l 5 



(12) 

If A„ o is the largest positive eigenvalue of W, ea. (|12[l then 
has a solution £ p = u>„ (p) <2i, , with 
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,(k)| 2 |Q,J 2 
(13) 

It is useful to compare ea. lfH3|) with eq.@. Since £ k 
is gapless, the integral in eq.© is logarithmically di- 
vergent as A — * 0. As a result, eq.® can always 
be satisfied (for all g) by an appropriate choice of A. 
Ea.ljT3")l is different. Due to the gap in E p , the integral 
in eq. I|13[) at h = is bounded by the convergent integral 
I Vo = tt- 1 J2 P \w Vo (p)\ 2 /E p . We then have the following 
possibilities : 

(i) For \ Vo > A c = l/I Vo , eq.JTSl implies Q Vo = C + 
Dh 2 + .. as h — * 0, where C and D are constants. In 
other words, for sufficiently strong interactions, triplet 
fluctuations (C) can exist in the the ground state of equal 
spin population. Ea. (|ll|l then implies to = M/Cl = + 
0(h 2 ), where \ is the susceptibility 



X(M,A) 



-y 
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and — W i/„0 S -)Q° i s the solution of ea. H13|) at h = 0. 

(ii) For X Vo < A c , Q Vo vanishes at h = 0. However, 
since the integral in ea. (|13fl is logarithmicaly divergent 
as h — ► A and Q„ o — ► 0, a nonzero Q Vo can always be 
found by increasing h before it reaches A. The critical 
h (denoted as h c ) when Q v becomes non-zero is given 
by A" 1 = n- 1 E k>CT =± |^ (k)| 2 /(£ k + ah c ). As h-,h c 
from above, we have <2„ oc (h— h^ 1 / 2 , m{h) oc (h — h c ). 

(iii) At A = A c , we have Q Vg ~ h and to ~ ft 3 . 

At unitarity, where the only energy scale at T — 
is /1 (or Fermi energy Ep), the system acquires a uni- 
versal thermodynamics [llj with magnetization given by 
m{jj) = n(fj,)G(h/ /1), where cx /1 3 / 2 is the number 
density, and G(h//j,) is a dimcnsionless universal func- 
tion. There is no need to consider the A-dependence 
of to, since at unitarity we have A = a/i, where a is a 



universal constant, (a — 1.16 in BCS theory, and 1.22 ac- 
cording to Quantum Monte Carlo calculations^^)- F° r 
large quasiparticle interactions, A > A c , to is linear ft. 
This implies G(h/ (J,) oc ft//i, or m(/x) oc y}l 2 h. In the 
presence of a trap V(r) (which vanishes at, say, r = 0), 
local density approximation (p,—*/j,— V(r)) implies 



m(r) = m(0)(l - V(r)/fi) 1/2 



(15) 



in the region where /i — V(r) > 0, a property that can 
be tested experimentally. For weak quasi-particle inter- 
actions so that G(h/n) is not linear in h for small ft, the 
full scaling form m((i) = n(fJ,)G(h/ '/z) has to be used to 
determine the spin density profile. 

C. A specific model of quasi-particle interaction 
at unitarity: In principle, the quasi-particle interac- 
tion is determined by the vertex function r(pi,p2;P3,P4)- 
As in studying fcrmion pairing, one needs to first 
identify the relevant interactions. Since spin fluctua- 
tions are equivalent to quasi-particle pairing (see Sec. 
A), we look for terms in the interactions that gen- 
erate triplet pairs with zero total momentum in the 
particle- hole channel[l3j, which are of the form V qp = 

fi_1 Ep,p'- F (P I P') c P 'T c Pl c P'l c PT' where F (P'P') an in " 
teraction function symmetric in p and p'. Expressing the 

c's in terms of a k and /3k, we find V qp = V qp + U qp where 

V qp is precisely eq.© with W(p,p') = F(p,p')x p X p >, 

with Xp = u p v P ■ The U qp are terms that do not produce 

triplet fluctuations and will be ignored|14j. 

Still, there is the problem of finding the expression of 

^Xp>pO- Once again, simplifications can be obtain at 

unitarity, since the only momentum scale at T = is 



k Q = \2M^/h 2 . If we further assume that the existence 
of an effective field theory for the fermions where the 
interaction coefficients can be expanded in powers of p, 
then in the lowest order in momentum, we have 



W(p,p') = F(p,p') Xp x P - 



4nh 2 Xp ■ p' 



Mk s -X P X pl - (16) 
where A is dimensionless. Ea. (|12fl then has a solution 
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FIG. 1: The triplet fluctuation A (in ea. <ll7t ') as a function of 
h/ [i. For A > A c , A = A c , and A < A c , we have A oc constant, 
A oc h and A oc ^Jh — h c respectively. 
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FIG. 2: The polarization m/n as a function of h/fi according 
to eg. IIH . For A > A c , A = A c , and A < A c , we have m oc h, 
m oc h z and m oc (ft — h c ) respectively. 
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FIG. 3: The difference of angular averaged momentum dis- 
tribution between different spin populations, Sn(k). 



C k = MXkP ' A / fc o> witn 

A= 4^1_ ^ (X k k)(k-A) 

Mk o ^ k ^ ± 2^ k -^) 2 + lx k k-A|2 

As discussed in Sec.B2, the solution of ea. l|17|) depends 
on whether A > A c , where X~ l = E k (x£l f k ' A) . 

o k 

Taking the BCS result A = 1.16/x at unitarity, we have 
A c = 2.14. We also find that the solution of ea.lfTTjl with 
lowest energy is A = ^4(x ± iy)/y/2. 

The behavior of triplet fluctuation A and polarization 
m/n as a function of h/\x for different A are shown in 
FiglD and [3 which illustrate the behaviors discussed in 
(i) to (iii) mentioned in Sec.B3. (Note that A oc Q v .) 
FigEl shows the difference between momentum distri- 
butions of different spin components, averaged over all 
angles, Sn(k) = (in)^ 1 J dk[n^(k) — n.j(k)] for differ- 
ent quasiparticle interaction A and different polarization 
m/n. It is interesting to note that the difference expands 
over a wide range of momenta around k a . Moreover, 
the difference is not very sensitive to whether A > A c or 
A < A c . For m/n = 0.1, the difference at k can be as 
large are 0.14. Measurement of momentum distribution 
can therefore be used to detect the existence of this new 
state eq.ljSJl. In addition to the profile for magnetiza- 
tion (eq.@5j)) and 5n{k), the new state eq.© can also be 
detected by noise measurements. Since eq.© contains 
It and || pairs, detection of correlations between like 
spin fermions with opposite momenta at small polariza- 
tion will make a strong case for the new state eq.©. 

Concluding remarks: We have shown that quasi- 
particle interactions can introduce triplet fluctuations in 
the BCS state. The resulting state eq.ljEJl will have non- 
zero "spin susceptibility" even at very low "magnetic 
fields" h. Although our calculations are based on a spe- 
cific model on the quasiparticle interaction, many con- 
sequences of these interactions can be deduced on gen- 
eral grounds. To compare energies with phase separation 
state requires a better understanding of the vertex func- 
tion T(pi,p2;p3,P4:) at unitarity, which is an intriguing 
and challenging problem in itself and will be studied else- 



where. In any case, the new state eq.JHJ) has many unique 
features in the density profile (eci. (|15|) ). in the momen- 
tum distribution (figure 3) , and noise correlation between 
like spins. All these properties can be measured experi- 
mentally. 
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